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ABSTRACT 

Due  to  quantum  confinement,  superlattice  detectors 
may  conceivably  be  tuned  to  detect  electromagnetic 
signatures  difficult  to  detect  with  technology  available 
today.  The  challenge  in  creating  new  superlattice 
detectors  is  in  the  growth  process  where  the  confluence  of 
chemistry,  physics,  and  mechanics  make  control  of 
parameters  nontrivial.  Furthermore,  the  time  scales 
involved  are  often  prohibitive  for  deterministic  modeling 
approaches.  We  have  developed  and  validated  a  new 
modeling  method  for  probabilistic  modeling  of 
superlattice  growth  that  spans  the  appropriate  time  scales 
relevant  for  experiments.  The  predictions  of  the  models 
have  been  tested  successfully  against  numerous 
independent  experiments.  Through  a  fundamentally  new 
Green's  function  formulation  for  strain  interactions 
among  quantum  dots,  we  present  a  new  kinetic  Monte 
Carlo  methodology  that  can  successfully  predict  pattern 
formation  on  surfaces  while  requiring  only 
experimentally-viable  parameters  as  input.  The  new 
method  also  can  be  used  for  engineering  design  of  dot 
correlations  in-plane  and  in-bulk  to  predict  dot  alignment 
in  general  superlattices. 

INTRODUCTION 

Devices  such  as  Fourier  Transform  Infrared  (FTIR) 
spectrometers  can  detect  electromagnetic  signatures 
related  to  the  energies  of  bonds  in  organic  molecules. 
These  energies  are  often  in  the  range  of  400  -  5000  cm"1 
representing  the  typical  absorption  spectroscopy  peaks  for 
specific  types  of  bonds.  For  example,  the  C-Br  bond  in 
bromoalkanes  has  an  accompanying  peak  at  500  cm-1 
while  the  O-H  bond  in  most  alcohols  and  phenols  are 
around  3000  cm-1.  The  central  detector  element  is 
typically  a  collection  of  charge-coupled  devices  (CCD) 
that  is  each  comprised  of  a  thin  semiconductor  material 
layer  whose  electronic  band  gap  energy  uniquely 
determines  the  energy  range  of  electromagnetic  signal  it 
can  detect.  Modem  fabrication  techniques  involve  growth 
of  these  thin  layers  on  inexpensive  substrates  such  as 
silicon.  The  grown  layers  can  often  be  as  thin  as  a  few 
monolayers  in  height  but  are  often  no  larger  than  a 
micron.  At  these  small  length  scales,  the  interaction 


between  the  growth  and  the  substrate  is  nontrivial  to  the 
extent  that  the  mismatch  in  the  theoretical  lattice 
constants  can  appreciably  distort  the  overall  band  energy 
of  the  final  system  and  be  used  as  a  design  consideration 
for  targeted  device  performance.  Thus  it  is  not 
uncommon  to  have  a  material  grown  on  a  substrate  detect 
signals  that  the  free  standing  material  cannot. 

Fabrication  of  these  types  of  detector  materials 
involves  complex  atomistic  processes.  It  is  now  well- 
known  that  during  fabrication,  the  interactions  among  the 
atoms  in  the  growing  layer  lead  to  stark  features  on  the 
surface  known  as  self-assembled  quantum  dots  (QDs). 
Self-assembled  QDs  have  been  intensely  investigated  due 
to  observed  optical  and  electronic  properties  with 
potential  applications  in  a  broad  range  of  new 
optoelectronic  and  semiconductor  devices  (BresslerHill, 
e  al.,  1995;  Bimberg,  et  al.,  1998).  The  indium- 
arsenide/gallium-arsenide  (InAs/GaAs)  heterostmcture  is 
a  typical  example,  which  is  characterized  by  a  large  lattice 
mismatch  between  alternating  In  As  and  GaAs  layers.  At 
the  growth  surface,  the  material  undergoes  a  transition 
from  a  two-dimensional  (2D)  cluster  to  three-dimensional 
(3D)  islands  (Stangl  et  al.,  2004;  Joyce  and  Vvedensky, 
2004). 

Experimental  data  indicate  that  InAs/GaAs  island 
array  is  governed  by  thermodynamics  (Schukin  and 
Bimberg,  1999;  Pcheljakov  et  al.,  1997).  The  physics 
behind  the  transition  from  2D  cluster  to  3D  islands  has 
been  studied  recently.  For  example,  using  atomic  force 
microscopy,  (Arciprete  et  al.,  2006)  studied  how  kinetics 
drives  the  2D  to  3D  transition  in  InAs/GaAs(001) 
heterostructure.  The  transition  from  thermodynamically  to 
kinetically  controlled  QD  self-assembly  was  also  studied 
in  (Musikhin  et  al.,  2005),  both  experimentally  and 
theoretically. 

It  is  argued  that  strains  on  the  surface  change 
diffusivity  and  surface  mobility  (Zandvliet  and  Poelsema, 
1999;  Sage  et  al.,  2000;  Zoethout  et  al.,  2000;  Pao  and 
Srolovitz,  2006)  and  therefore  may  be  the  reason  for  the 
enhanced  evaporation  (Sun  et  al.,  2000).  The  epitaxial 
system  can  lower  the  free  energy  by  transferring  atoms 
from  the  island  edge  to  the  upper  layer,  because  the 


transition  leads  to  a  decrease  in  the  contact  area  between 
the  substrate  and  a  new  layer  (Arciprete  et  ah,  2006).  As 
such  atomic  transitions  to  the  upper  layers  lead  to  the 
strain  field  relaxation.  The  tradeoff  between  the  cost  of 
additional  surface  energy  and  the  gain  of  energy  due  to 
elastic  relaxation  is  the  very  driving  force  for  the 
transition  from  2D  cluster  to  3D  islands  (Kern  and  Muller, 
1997;  Chaparro  et  al.,  2000).  These  arguments  provide 
reason  enough  to  consider  situations  when  the  atom  hop 
probability  to  an  upper  layer  can  be  higher  than  that  to  a 
lower  layer.  The  kinetics  of  such  a  process  could  be 
described  by  adjusting  the  probability  for  atoms  to  hop  up 
or  hop  down,  which  could  subsequently  lead  to  2D  cluster 
or  3D  islands  self-assembly  (Arciprete  et  al.,  2006). 


(Lam  and  Vlachos,  2001;  Lung  et  al.,  2005),  an  up/down 
ratio  is  introduced. 


First,  the  2D  hopping  probability  of  an  atom  from 
one  lattice  site  to  a  nearest  or  next  nearest  neighbor  site  in 
the  (x,y)-plane  is  still  governed  by  the  Arrhenius  law 
enhanced  by  the  long-range  strain  energy  field  (Pan  et  al., 
2004;  Pan  et  al,  2006;  Nurminen  et  al.,  2001;  Larsson, 
2001).  The  Arrhenius  law  governing  the  hopping 
probability  is  give  by 
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While  some  simple  computational  approaches  for 
2D  to  3D  self-assembly  were  discussed  before  (i.e., 
(Brunev  et  al.,  2001)),  no  complete  3D  QD  self-assembly 
model  has  been  developed  in  which  the  growth  from  2D 
cluster  to  3D  islands  and  the  corresponding  island  size 
equalization  can  be  clearly  illustrated  during  the  growth 
process.  In  this  work  we  develop  a  fast  multiscale  3D 
kinetic  Monte  Carlo  (KMC)  QD  growth  model.  It  is 
developed  from  our  original  (x,y)-plane  growth  model 
(Pan  et  al.,  2004;  Pan  et  al.,  2006)  by  a  fast  algorithm  for 
the  long-range  strain  energy  calculation  and  by 
introducing  the  up/down  ratio  for  lateral  self-organization 
(Pan  et  al.,  2007). 


where  v0  is  the  attempt  frequency  (=1013s_1),  T  the 
temperature,  and  kB  the  Boltzmann’s  constant.  Also  in 
equation  (1),  Es  and  En  are  the  binding  energies  to  the 
surface  and  to  the  neighboring  atoms,  respectively. 
Finally,  Estr(x,y),  as  a  function  of  the  plane  coordinates 
(x,y),  is  the  energy  correction  from  the  long-range  strain 
field  due  to  the  lattice  mismatch  between  the  substrate 
and  the  deposited  material.  As  this  long-range  strain 
energy  needs  to  be  calculated  repeatedly  during  the 
simulation,  we  have  developed  a  fast  algorithm  by  pre¬ 
calculating  the  energy  along  a  unit  circle  and  interpolating 
the  energy  at  any  location  afterwards  (Pan  et  al,  2007;  Pan 
and  Yang,  2003). 


With  the  proposed  program,  we  have  successfully 
simulated  the  transition  of  the  2D  cluster  to  3D  island 
process.  Furthermore,  depositions  with  different  flux  rates 
are  studied  to  show  how  the  flux  rate  affects  the  island 
equalization  during  the  deposition  process. 

Our  model  has  also  shown  clearly  the  importance  of 
the  interruption  time  during  the  growth  of  QD  islands.  For 
the  sake  of  extending  these  models  to  applications 
requiring  strict  control  on  the  ordering,  we  also 
investigate  the  feasibility  of  using  interrupted  growth  to 
create  desired  patterns  on  surfaces.  This  obviates  the 
prevailing  practice  of  hand-machining  substrate 
morphology  or  seeding  the  surface  through  an  extra  step 
to  encourage  growth  at  specified  locations.  We  believe 
that  this  model  will  provide  an  attractive  means  for 
producing  predictably  ordered  nanostructures. 


We  calculate  the  binding  energy  to  the  neighboring 
atoms  in  the  (x,y)-plane,  En,  using  the  following 
algorithm:  We  assume  that  the  strength  of  a  single  nearest 
neighbor  bond  is  Eh  (=0.3eV),  and  it  is  reduced  by  a 

factor  of  a  (=l/  V2  )  for  the  next  nearest  neighbors.  To 
evaluate  the  diffusion  barrier,  the  binding  energy  at  the 
site  P0,  where  the  diffusing  atom  is  located,  is  calculated 
to  be 

Ep  =  nEb  +  amEb  (2) 

with  n<4  and  m<4  being,  respectively,  the  number  of 
nearest  and  next  nearest  atoms.  Similarly,  for  the  site  P\ 
where  the  atom  is  going  to  hop  to,  we  have 

Et  =g(n'Eh+am'Eb)  (3) 


MODEL  DESCRIPTION 

The  3D  layer-by-layer  KMC  growth  model  is 
developed  from  our  2D  (x,y)-plane  growth  model  (Pan  et 
al.,  2004;  Pan  et  al.,  2006).  As  in  our  previous  model,  the 
important  contribution  from  the  long-range  elastic  strain 
energy  is  included  using  a  fast  algorithm  based  on  the 
Green’s  function  method  (Pan,  2002).  Furthermore,  to 
account  for  the  out-of-plane  (up  and  down)  movement  of 
the  atoms,  which  is  also  called  desorption  and  adsorption 


where  nf< 4  and  m’<4  are,  respectively,  the  number  of 
nearest  and  next  nearest  atoms  at  the  new  site  Pu  and  g 
(=0.2)  describes  the  coupling  between  the  adjacent  lattice 
sites  (Meixner  et  al.,  2003).  Therefore,  the  overall  binding 
energy  En  caused  by  neighbor  interactions  for  an  atom  to 
jump  from  site  P0  to  site  P\  is  given  by  the  difference  of 
the  binding  energies  at  the  corresponding  lattice  sites 

Ea=(n-g-n')Eb+(m-g-m')aEb  (4) 


Second,  the  binding  energy  to  the  surface  of  the 
(x,y)-plane,  i.e.  Es ,  is  assumed  to  be  constant  (iss=1.3eV) 
in  our  previous  2D  (x,y)-plane  self-assembly  model  (Pan 
et  ah,  2004;  Pan  et  ah,  2006).  However,  in  the  3D  case 
(figure  1),  the  binding  energy  to  the  surface  depends  on 
the  surface  geometry  of  both  the  initial  and  final  positions 
of  the  atoms.  Based  on  recent  molecular  dynamics  (MD) 


calculation  (Montalenti  and  Voter,  2001)  and  KMC 
simulation  experience  (Pan  et  al.,  2004;  Pan  et  al.,  2006; 
Meixner  et  al,  2003)  in  2D,  we  therefore  propose  the 
following  simple  equation  for  the  3D  surface  binding 
energy  calculation. 

E.  ={]--s)Es0+{p-g-p')aE!tt+(q-g-q')a-aEsl){5) 


Figure  1  Schematic  illustration  of  3D  QD  self-assembly  model.  An  atom  on  top  of  the  substrate  surface  (x,y)-plane  in  (a), 
the  relative  locations  of  the  atom  grid  on  the  (x,y)-plane  in  (b),  the  corresponding  binding  energy  Es  related  to  the  in-plane 
locations  in  (c). 


where  Es 0  is  the  binding  energy  for  the  atom  exactly  under 
the  selected  adatom  (figure  1),  p,  q  are  number  of  nearest 
and  next  nearest  atoms  in  original  position  (n’<4,  m'<4), 
p\  q'  are  the  number  of  nearest  and  next  nearest  atoms  in 
the  new  position  (n’<4,  m’<4).  The  other  two  parameters, 
g  and  a ,  are  used  to  scale  the  energy  contribution  from 
the  atoms  in  the  first  and  second  squares  (figure  la,  lb). 
Assuming  that  the  maximum  surface  binding  energy  to  be 
1.3eV  as  before  for  the  2D  growth  simulation  (Pan  et  al., 
2004;  Pan  et  al,  2006)  which  means  that  all  the  positions 
under  this  adatom  are  occupied  by  atoms,  we  can  find 
J£,s0=0.28eV  from  equation  (5)  by  back-calculation. 

Finally,  as  discussed  in  the  introduction,  the  growth 
system  always  tries  to  decrease  its  free  energy  by  moving 
atoms  at  the  edge  to  upper  layer  to  form  3D  islands.  The 
possibility  for  an  atom  to  hop  to  a  higher  or  lower  layer 
depends  on  the  material  properties  and  growth  conditions 
(Arciprete  et  al.,  2006).  To  account  for  the  up  and  down 
jump  activities  of  the  atoms,  the  parameter  up/down  ratio 
p  is  introduced,  which  is  defined  as  the  ratio  of  the 
possibility  for  edge  atoms  hopping  up  to  that  hopping 
down.  In  other  words,  if  we  let  Puv  be  the  jump  up 
possibility  and  Tdown  the  jump  down  possibility,  then  the 
up/down  ratio  p  =JPup/JPdown.  It  is  further  remarked  that  the 
up/down  ratio  actually  reflects  the  balance  between 
surface  energy  increase  and  strain  energy  decrease.  This 
physical  activity  is  illustrated  in  figure  2  where  layer  1  is 
the  substrate  and  the  atom  “A”  is  on  the  edge  of  layer  2  of 
the  deposited  atoms  (i.e.,  within  the  QD  island).  It  is  also 
understandable  that  the  up  or  down  jump  possibility  of 
atom  “A”  is  controlled  only  by  the  up/down  ratio  instead 
of  by  the  individual  jump  up  and  down  possibilities. 
Moreover,  strictly  speaking,  since  the  strain  energy 
changes  along  the  island  height  direction,  the  up/down 
ratio  is,  in  general,  not  constant  in  different  island  layers. 


However,  in  order  to  extract  the  important  contribution  of 
the  up/down  ratio,  we  assume  that  the  up/down  ratio,  is 
constant.  In  other  words,  it  is  only  affected  by  the 
geometry  around  it  but  not  by  the  layer  position. 
Furthermore,  in  order  to  form  3D  islands,  the  number  of 
atoms  jumping  up  should  be  larger  than  those  jumping 
down,  which  means  that  p  should  be  larger  than  1. 
Otherwise,  all  the  atoms  will  tend  to  move  to  lower  layer 
to  form  the  layer-by-layer  Frank-van  der  Merwe  structure. 


Figure  2.  Illustration  of  the  up  or  down  jump  for  an  edge 
atom  A  during  3D  QD  self-assembly. 

RESULTS:  PARAMETER  STUDIES  OF  p 

Using  the  3D  QD  self-assembly  approach  described 
above,  we  can  now  simulate  the  epitaxial  growth  process. 
The  growth  model  contains  100x100  grids  with  periodic 
boundary  conditions,  and  the  material  is  InAs/GaAs 
(001).  We  first  study  the  effect  of  variations  in  p  on  the 
island  height.  The  validation  of  recognizable  pattern 
formations  was  performed  extensively  with  independent 
results  in  the  literature  in  (Pan  et  al.,  2007)  and  (Zhu  et 
al.,  2007).  Shown  in  figure  3  are  the  islands  distributions. 
It  is  observed  clearly  from  figure  3  that  1)  the  average 
island  height  increases  with  increasing  p  (for  small  p ,  say, 
p  <5,  basically  only  2D  growth  is  observed);  2)  island  size 
and  shape  are  very  sensitive  p  when  it  is  less  than  20;  3)  it 
is  worthy  to  remark  that  p  determines  the  2D  to  3D 


transition  instead  of  the  absolute  value  of  Pup  and  Pdown. 
This  ratio  is  directly  proportional  to  the  ratio  of  the 
ascending  (desorption)  and  descending  (adsorption)  atoms 
on  the  surface  (Lam  and  Vlachos,  2001;  Russo  and 
Smereka,  2006). 

To  further  understand  the  effect  of  variations  in  p  on 
the  average  island  height  the  data  in  figure  3  are  analyzed 
along  with  more  simulated  results  for  large  up/down 
ratios.  The  relation  between  the  average  island  height  and 
p  is  presented  in  figure  4.  Figure  4  demonstrates  that  1) 
the  average  island  height  experiences  sharp  changes  when 
p  varies  from  3  to  around  20.  The  island  height  increases 
with  increasing  p\  2)  when  p  is  larger  than  20,  the  curve 
asymptotically  approaches  a  maximal  island  height  of 
approximately  20  grid  spacings.  This  means  that  the 
average  island  height  will  mostly  keep  at  a  constant 
maximum  value  when  the  up/down  ratio  is  large  (say, 
p> 30);  3)  an  identifiable  inflection  point,  approximately  at 
p=  13,  physically  corresponds  to  the  critical  transition 
point  from  2D  cluster  growth  to  3D  islands  growth.  This 
up/down  ratio  is  roughly  equivalent  to  the  surface  and 
bulk  energy  ratio  as  was  demonstrated  in  (Brunev  et  al., 
2001;  Kratzer  et  al.,  2006)  using  different  approaches. 


Figure  3  3D  island  distributions  for  different  up/down 
ratios  p  with  total  coverage  c=1.6ML.  Fixed  growth 
parameters  are  T=700K,  FMhOlML/s,  and  interruption 
time  £i=200s  (The  total  simulation  time  =  the  deposition 
time  (^=160s)  plus  the  interruption  time  (^=200s)). 


Figure  4.  Average  island  height  vs.  p  for  p=  3  to  50.  Fixed 
growth  parameters  are  T- 700K,  F=0.01ML/s,  c=1.6ML, 
and  interruption  time  ^=200s. 

RESULTS:  ISLAND  SIZE  &  DISTRIBUTION 

Island  size  equalization  relies  on  the  movement  of 
atoms  on  the  surface  during  epitaxial  growth. 
Equalization  starts  from  the  very  beginning  of  the 
deposition  process  and  reaches  equilibrium  with  sufficient 
growth  time  (Pan  et  al.,  2004).  Figure  5  demonstrates  the 
island  size  dependence  on  growth  times  for  three  different 
up/down  ratios  (p=  10,  20,  and  30  in  the  first,  second,  and 
third  row).  Island  distributions  immediately  after  the 
deposition  (i.e.,  the  interruption  time  /i=0)  are  shown  in 
the  first  column,  100  seconds  after  the  deposition 
(^i=100s)  in  the  second  column,  and  200  seconds  after 
deposition  (^=200s)  in  the  third  column  (figure  5).  It  is 
observed  from  figure  5  that:  1)  for  a  fixed  up/down  ratio 
(p=  10,  20,  and  30),  with  increasing  interruption  time, 
small  isolated  islands  assemble  to  form  large  ones,  and 
the  island  distribution  becomes  increasingly  ordered;  2) 
for  fixed  interruption  time  (^=0,  100s,  and  200s),  an 
increasing  up/down  ratio,  in  general,  increases  the  height 
of  the  islands  (changes  the  shape  of  the  islands),  similarly 
observed  in  figure  3.  Since  p  is  closely  related  to  the 
surface  and  bulk  energy  ratio,  this  may  enable 
experimentalists  with  additional  control  to  optimize  the 
island  shape  and  distribution. 

While  figure  5  illustrates  the  influence  of  the 
interruption  time  on  the  island  size,  figures  6  and  7 
demonstrate  further  the  effect  of  the  deposition  time  on 
the  island  shape  and  size  for  fixed  up/down  ratio  p=  10.  In 
figure  6,  island  size  equalization  starts  from  the  beginning 
of  the  deposition  process  and  becomes  larger  with 
increasing  deposition  (£d=50s,  100s,  160s)  and 

interruption  time  (^=200s).  Furthermore,  with  increasing 
island  size,  the  number  of  islands  decreases.  This  relation 
can  be  observed  from  the  histogram  in  figure  7  for  the 
relationship  between  the  island  diameter  (of  the 


equivalent  circle  at  the  bottom  of  the  island)  and  the 
corresponding  number  of  islands. 


Figure  8  shows  the  effect  of  the  flux  rate  (F=lMl/s, 
O.lML/s  and  O.OlML/s)  on  the  island  distribution.  For  a 
low  flux  rate  (e.g.  at  F=0.01Ml/s),  the  atoms  are  afforded 
more  time  to  move  to  the  equilibrium  position  during  the 
deposition  process  and  to  assemble  together.  A  low  flux 
rate  (say,  F=0.01  ML/s)  usually  corresponds  to  larger  and 
more  ordered  islands. 


Figure  5.  Island  distributions  for  different  interruption 
times  (£i=0s,  corresponding  to  deposition  time  £d=160s, 
£i=100s,  and  ^=200s)  and  for  different  up/down  ratios 
(p=  10,  20  and  30).  Fixed  growth  parameters  are  T=700K, 
F=0.01ML/s,  and  c=1.6ML. 
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Figure  6.  Island  distributions  during  and  after  deposition. 
Deposition  time  £d=50s  with  0.5ML  coverage  in  (a), 
£d=100s  with  1ML  coverage  in  (b),  £d=160s  with  coverage 
1.6ML  in  (c)  and  interruption  time  ^=200s  in  (d). 
T=700K,  F=0.01ML/s,  c=1.6ML,  and  p=  10. 


RESULTS:  PREPATTERNING 


(Atkinson  et  al.,  2006)  or  masked  substrates  (Liang  et  al., 
2004),  both  of  which  require  substantial  increase  in  total 
fabrication  time.  A  nominal  change  to  the  interruption 
time,  however,  reveals  that  pattern  quality  (i.e.  larger  size 
and  repeated  ordering  of  dots)  can  be  obtained  (Pan  et  al., 
2007).  In  figure  9,  significant  differences  in  pattern 
quality  can  be  observed  between  traditional  continuous 
deposition  (figure  9a)  versus  employing  a  growth  pause 
for  50  seconds  after  an  initial  0.3ML  deposit  and 
interruption  at  150  seconds  after  1.6ML  deposit. 


Figure  7.  Average  island  diameter  (of  the  bottom 
equivalent  circle  of  the  island)  vs.  number  of  islands 
during  and  after  deposition.  T=700K,  F=0.01ML/s, 
c=1.6ML  andp=10. 
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Figure  8.  Island  distributions  for  different  deposition 
processes  of  0.5ML,  1ML  andl.6  ML  with  flux  rates 
F=lMl/s,  O.lML/s  and  O.OlML/s.  Fixed  growth 
parameters  are  T=700K,  c=1.6ML  andp=10. 


Control  of  pattern  formations  is  currently  enabled 
through  premachining  of  the  substrate  with  pore  arrays 


RESULTS:  CORRELATED  GROWTH 


Figure  9.  Comparison  of  (a)  continuous  growth  and  (b) 
prepatterned  growth  using  variations  of  substrate 
symmetries  (isotropic,  (001),  and  (111))  at  T=700K, 
coverage  of  1.6ML  and  p=10. 


Finally,  using  the  fundamental  developments  from 
(Pan,  2002)  that  are  used  in  the  3D  KMC  model,  the 
effect  of  the  buried  height  of  nanostructures  on  correlated 
growth  in  superlattices  is  investigated.  It  has  been  shown 
that  the  depth  to  which  QDs  are  buried  in  layered 
structures  influences  patterns  on  subsequently  grown 
layers  (Holy  et  al.,  1999;  Wang  et  ah,  2004).  The 
mechanism  that  enables  this  is  the  interaction  of 
mechanical  strain  fields  owing  to  the  embedding  strain  of 
the  array  of  QDs.  In  particular,  Wang  et  al.  showed  the 
strain  driven  correlation  and  anticorrelation  of  InGaAs 
QDs  in  detail  through  cross-sectional  STM  images. 

Using  the  Green’s  Functions  developed  in  (Pan, 
2002),  we  computed  the  strain  energy  values  on  the 
surface  due  to  a  buried  structure  in  an  elastic  half-space. 
The  degree  of  burial  is  varied  to  modulate  the  degree  of 
strain  interaction  between  the  buried  QDs  and  the  surface. 
The  dark  regions  in  Figure  10  show  regions  of  increased 
strain  energy  where  the  higher  accompanying  strain 
implies  that  coherent  structures  like  QDs  are  more  likely 
to  grow.  This  is  due  to  the  local  relaxation  of  the  degree 
of  lattice  mismatch  at  those  regions  where  the  strain  fields 
interact  constructively.  It  is  noteworthy  that  for  D=2  the 
(001)  growth  is  correlated,  which  agrees  convincingly 
with  those  of  the  STM  images  in  (Wang  et  al.,  2004)  at 
small  burial  heights.  At  D=20  the  (001)  growth  is  anti¬ 
correlated  which  also  corresponds  to  the  anti-correlation 
achieved  with  increased  burial  depth  in  (Wang,  et  al., 
2004). 

It  is  furthermore  noteworthy  that  from  the 
calculations,  we  show  that  the  (111)  growths  do  not 
exhibit  anti-correlation  by  increasing  the  burial  depth  to 
D=20  in  contrast  with  the  observations  from  the  (001) 
case.  In  fact,  the  growth  remains  strictly  correlated. 
Although  the  (111)  growth  correlation  remains  to  be 
shown  experimentally,  the  absence  of  anti-correlation  at 
deeper  burial  depths  appears  to  be  caused  by  the  lack  of 
strain  field  interactions  at  the  anti-correlated  sites  due  to 
the  additional  effect  of  anisotropy  and  the  overall  change 
in  symmetry  caused  by  the  change  in  growth  direction. 


Figure  10.  Strain  energy  contours  of  3x3  buried  array  at  two 
areas  indicate  regions  of  favorable  growth  for  the  next  layer. 

CONCLUSIONS 

In  this  report,  we  propose  a  3D  QD  epitaxial  growth 
model  by  enhancing  our  former  (x,y)-plane  growth  model 
with  a  fast  algorithm  for  long-range  strain  energy 
calculation  and  an  up/down  jump  ratio  for  the  diffusion 
atoms.  Specifically,  the  balance  between  surface  energy 
increase  and  bulk  energy  decrease  is  demonstrated  by 
introducing  the  up/down  ratio.  Further  studied  is  the 
dependence  of  the  island  height  and  shape  on  the  up/down 
ratio.  Combining  the  up/down  ratio  and  one  of  the  growth 
parameters  (i.e.,  flux  rate),  island  equalization  during 
deposition  process  and  after  deposition  is  demonstrated. 

The  developments  have  also  enabled  three  new 
fundamental  results.  The  first  is  the  consistent  validation 
of  pattern  formations  for  general  quantum  dot  growth  on 
surfaces.  By  considering  only  the  chemical  species  used 
in  the  material  fabrication,  fabrication  parameters  such  as 
temperature  and  flow  rate,  and  the  symmetry  of  the 
substrate  surface  plane,  we  can  repeatedly  reproduce  the 
patterns  obtained  from  numerous  independent  results 
from  the  literature.  Based  only  on  the  experimental 
parameters  reported  in  those  papers,  the  new  modeling 
theory  accurately  reproduces  the  quantum  dot 
arrangements. 


different  depths  and  two  different  growth  directions.  Dark 


The  second  result  is  a  new  fabrication  procedure  that 
obviates  tooling  for  prepatterned  growth  when  strict 
control  is  required  on  the  first  surface  grown  directly  on 
the  substrate  layer.  By  adjusting  species  flow  rates  during 
growth,  we  have  simulated  a  technique  for  controlling  the 
pattern  quality  on  surfaces.  Whereas  earlier  efforts 
attempted  to  seed  growth  by  mechanically  creating 
“dimples”  on  the  substrate,  the  new  method  enables 
ripening  to  coalesce  smaller  dots  into  larger  ones,  thereby 
producing  more  pronounced  dot  sizes  in  specifiable 
patterns  without  machining  the  substrate. 

The  third  result  is  the  development  of  a  through¬ 
thickness  dot  correlation  technique  for  superlattices  where 
elastic  fields  generated  from  buried  quantum  dots  are 
demonstrated  to  provide  another  means  of  controlling 
surface  pattern  quality.  Dots  among  different  layers  can 
interact  and  produce  correlated  and  anticorrelated 
patterns.  Although  the  capping  layer  thickness  has  been 
known  to  produce  this  effect,  these  simulated  results  are 
the  first  of  its  kind. 
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